Mechanism for Powerlaws without Self-Organization 
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Abstract: A recent claim has been made that there must be a self-regulation in the waiting times to 
see hospital consultants on the ground that the relative changes in the size of waiting lists follow a power 
law [Q]. In agreement with simulations of Frecketon and Sutherland, we explain the general non-self- 
regulating mechanism underlying this result and derive the exponent value —2 exactly. In addition, we 
provide links with related phenomena encountered in many other fields. 
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Power law probability distribution (pdf) functions of event or object sizes seem to be an ubiquitous 
statistical features of natural and social systems. It has repeatedly argued that such an observation calls 
for an underlying self-organizing mechanism. In the last decade, such claims have been made for instance 
for earthquakes, weather and climate changes, solar flares, the fossil record, and many other systems, to 
invoke the relevance of self-organized criticality [|l|]. This claim is often unwarranted as there are many 
non-self-organizing mechanisms producing power law distributions ||^ ^]. Recently, a similar claim of 
the existence of self-regulation in the waiting time to see hospital consultants has been made on the 
ground that the relative changes in the size of waiting lists follow a power law [Q]. Using numerical 
simulations, Frecketon and Sutherland ^ show however that such a power law for relative changes can 
be reproduced by the simple non-self-regulating random walk null-hypothesis, which obtains a power 
law exponent (slope of the log-log plot of the non-cumulative distribution) indistinguishable from —2. 

Here, we explain the general non-self-regulating mechanism underlying this result and derive the 
exponent value —2 exactly. This mechanism is nothing but the effect of the conservation of probability 
Px{x) dx = Pyiy) dy on the probability distribution function (pdf) Px{x) of a change of a variable x 
to another y = f{x). In other words, the objective estimation of the probability of an event is invariant 
under a change of mathematical description of this event. If the transformation is an inversion y = x^^^" 
with a > 0, then 
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Suppose that Px{x) goes to a constant for x ^ 0, then the distribution of y for large y is a power law 
with a tail exponent 1 + a. The uniform fluctuations of x close to zero lead to scale-free and arbitrarily 
large fluctuations of its inverse power y. The power law form is kept obviously (with a modification of 
the value of the exponent) if Px{x) itself goes to zero or diverges close to the origin as a power law. 

This mechanism operates in many situations. Let us illustrate it in details in a simple well-controlled 
situation to demonstrate that it also applies to rather complex systems. Consider the Ising model of 
interacting spins with ferromagnetic interactions which tend to align the spins while thermal fluctuations 
tend to randomize them. The Ising model is a paradigm of the fight between order and disorder leading 
to a spontaneous organization by collective behavior and has been used to describe both physical and 
social phenomena. Here, we recall that the distribution of the relative changes AM/M of magnetization 
M (total normalized sum of the spin values) is a power law pdf, but this power law owes nothing to 
organization and everything to the "change-of-variable" mechanism with an inversion. Right at the 
critical temperature T = Tc, the pdf P{M) of magnetization has a very "thin" exponential tail P{M) oc 
—const M^+^}, with 5 + 1 « 5.8 for the 3D Ising universality class [^]. The distribution 
of changes AM of the magnetization M under a fixed number of Monte Carlo steps per site is a Gaussian 
of standard deviation 1^ at T = Tc ^J. However, the distribution of the relative changes AM/M is a 
power law with exponent — 2 [^]! How does it come about? 

Large values of AM/M come from the limit M rather than AM —>■ oo. The probability P{M) 
goes to a constant for M ^ 0, while AM can be approximated by the width W of the Gaussian; thus: 
P{X)dX = P{AM/M)d{AM/M) = P{M)dM ^ const dM and P{X) = P{M)/{dX/dM) = 
const /{d{W/M)/dM) oc in agreement with Jan et al. |^]. The l/X"^ power law is in fact not 

restricted to the critical point and is very general since it results simply from the inversion mechanism 
y = x"^/" defined above with a = 1. A similar argument has been proposed in in this context. 

The power law pdf for the random walk null-hypothesis studied by Frecketon and Sutherland ||5|] 
results from the choice of the variable which is the relative change in numbers {Nt+i — Nt)/Nt. Ac- 
cording to the "change-of-variable" mechanism with an inversion, the power law tail is controlled by 
excursions of the denominator towards small values. Since again a = 1, this predicts the value —2 for 
the exponent in agreement with the numerical simulations It is tempting to conclude that the same 
general mechanism operates for the results of Smethurst and Williams [Q]. 

The same "change-of-variable" mechanism with an inversion applies in many other situations to 
create power law pdf 's (see for details and references): 

• In continuous percolation, the pdf of transport coefficients such as conductance, permeability and 



rupture thresholds and of necks between random holes or random cracks [10, 11, N 



The pdf of fluid velocities due to the presence of vortices in hydrodynamics |[T^ ; 

The Holtsmark's distribution of gravitional forces created by a random distribution of stars in an 
infinite universe (see [14] and chapter 17 of [||] and references therein); 

This result for the Holtsmark's distribution applies to other fields with a suitable modification of 
the exponent, such as electric, elastic or hydrodynamics, with a singular power law dependence of 
the force as a function of the distance to the source (see chapter 17 of [^] and references therein). 

The statistics 

[(n-l)-iE-=i(x,-x)2]i/2 ^ 

occurs in the construction of tests and confidence intervals of the random variable and is described 
by Student's distribution, whose power law tail results from the same "change-of-variable" mech- 
anism with an inversion. Here the inversion refers to the fact that the estimated standard deviation 



in the denominator of T can approach zero leading again to large fluctuations of T. For a demon- 
stration, see chapter 14 of [^] and references therein. 



Berry's "battles of catastrophes" [15]; 



-7/2 power law pdf of density in the Burgers/adhesion model of the universe [16]. 



In sum, this communication has explained one among many other mechanisms for power law pdf 's 
(see chapter 14 of ^ for a general review) that are not based on self-regulation or self-organization. By 
this understanding of the main possible mechanisms, and complementing this toolbox with mechanisms 
relying on collective effects, such as in percolation, criticality and self-organized criticality, it should 
become possible to identify which one is the most relevant to a given problem. 
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